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SUMMARY 

The Institute for Computational Mechanics in Propulsion (ICOMP) was formed to develop techniques to improve problem-solving 
capabilities in all aspects of computational mechanics related to propulsion. ICOMP is operated by the Ohio Aerospace Institute (OAI) 
and funded via numerous cooperative agreements by the NASA Lewis Research Center in Cleveland, Ohio. This report describes the 
activities at ICOMP during 1998, the Institutes thirteenth year of operation. 


INTRODUCTION 

The Institute for Computational Mechanics in Propulsion (ICOMP) was established at the NASA Lewis Research Center in September 
1985. The overall purpose was to improve problem-solving capabilities in all aspects of computational mechanics relating to propulsion. 
ICOMP provides a means for researchers with experience and expertise to spend time in residence at Lewis performing research to 
improve computational capability in the many broad and interacting disciplines of interest in aerospace propulsion. 

The scope of the ICOMP program is to advance the understanding of aerospace propulsion physical phenomena and to improve computer 
simulation of aerospace propulsion systems and components. The specific areas of interest in computational research include, fluid 
mechanics for internal flows; CFD methods; turbulence modelling; and computational aeroacoustics. 

This report summarizes the activities at ICOMP during 1998. 

The following sections of this report provide lists of the resident and visiting researchers, their affiliations and educational backgrounds. 
Individual sections are provided which briefly describe reports of RESEARCH IN PROGRESS, the REPORTS AND ABSTRACTS 
published over the past year. 
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K QM P R ES EAR C H S TA FF 


THE ICOMP STAFF OF VISITING RESEARCHERS 

The ICOMP research staff for 1998 is shown in Table I. A total of eighteen researchers were in residence at ICOMP for periods varying 
from a few days to a year. The resident staff numbered fifteen while the visiting staff numbered three. 

As usual, the resident researchers were very productive Table II provides a numerical summary of ICOMP during its thirteen years of 
operation in terms of research staff size and technical output as measured by the numbers of seminars, reports and workshops. 
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Table I. - The ICOMP Research Staff-1998 


A. Resident Staff. 

Joongkec Chungs Ph£>., Mechanical Engineering, University of California, Berkeley, 1991. May, 1992-Present. 

Duane R Hixon . PhD., Aerospace Engineering, Georgia Institute of Technology, 1993. October, 1993-Present. 

Tames M. Loellbach. Ph.D. expected 1995, Aeronautical and Astronautical Engineering, University of Illinois. May, 1992-Present. 
Andrew T. Norris. PhD., Mechanical and Aerospace Engineering, Cornell University, 1993. June, 1993-Present. 

Knsh Radhaknshnan. Ph.D., Mechanical Engineering, Massachusetts Institute of Technology, 1979. January, 1998-Present. 

A amir Shabbir. Ph.D., Mechanical Engineering, State University of New YorkJBuffalo, 1987. June, 1991-Present. 

Scott Sherer. M.S., Aeronautical and Astronautical Engineering, 1994. September, 1997-September, 1998. 

Shvue-Homg Shih. PhD., Aerospace Engineering, University of Cincinnati, 1993. June, 1993-October, 1998. 

Tsan-Hsing Shih . PhD., Aerospace Engineering, Cornell University, 1984. March, 1990-Present 

Erlendur Steinthorsson. PhD., Mechanical Engineering, Carnegie Mellon University, 1992. January, 1992-October, 1998. 

Gerald Tmmmer . B.S., Michigan State University, 1982. September, 1995-Present. 

Fu-Lin Tsung. Ph.D., Aerospace Engineering, Georgia Institute of Technology. March, 1993-Present. 

Michael White. PhD., Mechanical Engineering, University of California, Davis. October, 1997-Present 
Shave Yungs ter. PhD., Aeronautics and Astronautics, University of Washington, 1989. November, 1989-Present. 

Ge-Chcng Zha. PhD., Mechanical Engineering, University of Montreal, 1994. July, 1997-October, 1998. 

B. Visiting Staff/ Consultants. 

Thomas Haystrom. Ph-D., Applied Mathematics, California Institute of Technology, 1983. Associate Professor, Department of 
Mathematics and Statistics, University of New Mexico. 

Rli Turkel. Ph.D., Applied Mathematics, New York University, 1970. Professor, Department of Mathematics, Tel Aviv University, Tel 
Aviv, Israel. 

Mei Zhuang. Ph.D., Professor, Department of Mechanical Engineering, Michigan State University. 
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RESEARCH IN PROGRESS 


JOONGKEE CHUNG 

Research Anas: Cede Development for Unsteady Inlet Flows Using Parallel Processing, Iced Airfoils 

| ^ experiment was conducted in the Icing Research Tunnel at NASA Lewis to deter m ine the ice shapes which might accrete on the 

wings of a turboprop aircraft that had been in an icing accident This experiment was performed as part of NASA Lewis's support 
of the National Transportation Safety Board (NTSB) aircraft accident investigation. Using a vertically mounted 6 foot wing cut from 
the actual aircraft, v ario us ice shapes, depending on the ici n g conditions, were acreted on the wing. Subsequently, these shapes were 
to be used in an extensive numerical analysis. The experimental model had inflatable boots which could be activated during the 
experiment to shed the ice shapes. 

2. As a first step to investigate performance degradation due to icing through numerical investigation, a mathematical model using 
surface smoothing of the grid was applied. This approach was implemented in an interactive code, TURBO -GRD to generate surface 
shapes using different levels of smoothing. The code constructs a smooth curve whose shape is controlled by a piece-wise linear curve 
formed from selected discretized points. This study showed that a 50% level of CPS or higher is required to adequately represent 
the ice shapes as measured by having marginal influence, Le., less than 5% variation on the resulting lift and drag values. The details 
of the investigation were reported in ALAA Paper 98-3242 entided "Effect of Airfoil Ice Shape Smoothing on the Aerodynamic 
Performances” the paper was presented at the AIAA/ASME/SAE/ ASEE Joint Propulsion Conference held in Cleveland in July, 
1998. 

3. A numerical study was also performed to determine the aerodynamic effects of ice contamination on the wing of a commercial 
turbo-prop-powered aircraft Based upon the results of the study mentioned above, a prominent ndged-ice formation on the subject 
aircraft wing was selected for detailed flow analysis using 2-D, as well as, 3-D Navier-Stokes computations. This configuration was 
selected because it has die largest lift decrease and drag increase among all the ice shapes investigated in the study. A gnd sensitivity 
test was performed to determine the influence of grid spacing on the lift, drag, and associated angle-of-attack for the maximum lift. 
This study showed that grid resolution is important and a sensitivity analysis is an essential element of the process in order to assure 
that die final solution is grid independent The 2-D results suggested that a severe instability and an associated control difficulty could 
have occurred at a slighdy higher angle-of-attack than the one recorded by die Flight Data Recorder. This instability and control 
problem was thought to have resulted from decreased differential lift on the wings with respect to the normal loading for the 
configuration. The analysis also indicated that this instability and control problem could have occurred whether or not natural ice 
shedding took place. Numerical results using an assumed 3-D ice shape showed an increase of about 4 degrees in the angle at which 
this phenomena occurred when this occurred using a 2-D solver, trailing edge separation was also observed but started only when 
die angle of attack was very dose to the angle at which the maximum lift occurred. Detailed results of the experimental and numcncal 
study of this accident investigation were presented at the 37th ALAA Aerospace Sciences Meeting in January 1999. The tides and 
the authors of the papers are: 

(1) A. Reehorst, J. Chung, M. Potapczuk, Y. Choo, W. Wnght, and T. Langhals, "An Experimental and Numerical Study of Icing 
Effects on the Performance and Controllability of a Twin Engine Aircraft", AIAA Paper 99-0374, Reno, NV, Jan., 1999. 

(2) J. Chung, Y. Choo, A. Reehorst, M. Potapczuk, and J. Slater, "Navier-Stokes Analysis of the Flowfield Characteristics of an 
Ice Contaminated Aircraft Wing", AIAA Paper 99-0375, Reno, NV, Jan., 1999 

4 As preparation for an experimental investigation using a Natural-Laminar-Airfoil in the Low Turbulence Pressure Tunnel at NASA 
Langley, a series of numerical analysis were performed to hdp reduce the test matrix, to aid the maximum stress analysis, etc. The 
surface shapes of the stainless steel piece for the pressure measurement were generated using the mathematical model described 
above. It is estimated that half of the problem could have occurred whether or not natural ice shedding took place. 

THOMAS HAGSTROM 

Research Areas: Algorithms for Boundary Layer Value Problems , Domain Decomposition 
Research work during FY98 involved two efforts: 

Radiation Boundary Conditions for Computa tional Acroacousrics 

Recent years have seen a number of exciting developments in the practical implementation of high-order and exact radiation boundary 
conditions for the wave equation and Maxwell’s equations. In joint work with S. I. Hariharan of ICOMP we have developed an extremely 
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simple way to apply exact boundary conditions for die wave equation on a spherical boundary which uses only local operators. With other 
authors (B. Alpert and L Greengard), we have shown how to minimize the work in the application of nonlocal fo rmula tions, for both 
periodic/duct and exterior geometries, through the use of uniform rational approximants. Our ICOMP work has focused on the 
application of these and conditions developed earlier to die Euler equations linearized about a uniform flow. Particular accomplishments 
in 1997-98 include: 

L WithJ. Goodrich of die acoustics branch, a spatially periodic problem with an exact solution in integral form was developed and used 

to test the convergence of sequences of boundary conditions based on Pade approximants of increasing degree, and to compare the 
results with lower order treatments such as those proposed by Giles. As expected, for fixed times the Pade approximants exhibited 
‘spectral’ convergence, though the degree required for good accuracy increased linearly with increasing time. This behavior is 
consistent with the error estimates we have developed. It is worth noting that the problem chosen is quite challenging - results 
obtained with the Giles conditions were 40-80% in error after ten nondimensional time units. The test problem will also be used 
to examine the behavior of other boundary conditions such as stablized PMLs and unif orm rational approximants. We have also 
extended our exact solution to a duct geometry. 

i L We have extended die boundary conditions based on uniform rational approximants to the linearized Euler equations in periodic or 

duct geometry. These should prove more efficient than foe Pade conditions for long time computations, which may be relevant for 
cascade problems with time periodic forcing. 

iiL With S.I. H an h a r a n of ICOMP, we have extended our boundary conditions based on progressive wave expansions to the convective 
wave equation, preserving their property of ‘exactness’ on spherical boundaries for solutions described by finitely many spherical 
harmonics. It should be a simple step to go from the convective wave equation to the linearized Euler system. We have also made 
progress in generalizing the conditions based on uniform rational approximants (for exterior problems) to this case. 

iv. We have carried out preliminary work on die formulation and application of exact and approximate boundary conditions for the Euler 
equations linearized about parallel flows. 

High-Order Single Step Methods for Hyperbolic Systems 

J. Goodrich and R. Dyson of the acoustics branch have been developing and experimenting with high-order single step methods for solving 
hyperbolic systems. Particularly promising among these are methods based on Hcrmite interpolation. Experiments indicate they may be 
useful for simulations in complex geometries using Cartesian meshes, which is particularly efficient for grid generation and parallel 
computation. We have carried out preliminary work on the stability theory of the Hermite schemes, which we hope will provide further 
insight into the stability of various boundary treatments, and have also developed variable coefficient test problems to test various 
generalizations of the methods to the variable coefficient case. 

DUANE R. HDCON 

Research Area : Aeroacoustics 

Work continued on the development and application of high-accuracy numerical schemes for die computational prediction of supersonic 
jet noise. The sixth-order compact code has been rewritten in Fortran 90 to take advantage of the improved memory and data structures, 
and has been extensively validated The scheme and filters have been extensively tested on 1-D nonlinear benchmark problems and have 
shown very good results 11 - 21 . Tests were also conducted to assess the code's capability for complex geometries in preparation for a jet nozzle 
calculation 131 . A set of new boundary conditions were compared to see if any improvements could be made w . 

A collaborative effort has been underway with United Technologies Research Center to compare results from the Large Eddy Simulation 
codes, and initial results have been published 13,61 . These results compare the existing 2-4 MacCormack solver to UTRC's high-order upwind 
scheme; the new solver will be applied to this problem in 1999. 

Work will continue in two directions: application and improvement of the curvilinear code. The first 3-D LES nozzle + jet application 
will be run in 1999, and improvements in the accuracy, speed, and robustness of foe compact scheme will be investigated. The code will 
also be parallelized. 

References 

1) Hixon, R. (1998) Nonlinear Comparison of High-Order and Optimized Finite-Difference Schemes, ALAA Paper 98-3240, 34th 
AIAA/ASME/SAE Joint Propulsion Conference, Cleveland, OH. 
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2) Hixon, R. (1999) Prefactored Compact Filters for Computational Aeroacoustics, AIAA Paper 99-0358, 37th AIAA Aerospace 
Sciences Meeting, Reno, NV. 

3) Hixon, R, Shih, S. H., Mankbadi, R R, and Scott, J. R (1998) Time-Domain Solution of the Airfoil Gust Problem Using a High- 
OrderCompact Scheme, AIAA Paper 98-3241, 34th AIAA/ ASME/SAE Joint Propulsion Conference, Cleveland, OH. 

4) Hixon, R, Shih, S. H., and Mankbadi, R R (1998) Evaluation of Boundary Conditions for the Gust-Cascade Problem, AIAA Paper 
98-2489, 2nd AIAA Theoretical Fluid Mechanics Meeting, Albuquerque, NM Accepted for publication in Journal of Propulsion and 
Power. 

5) Shih. S. H., Hixon, R-, Povinelli, L. A., and Mankbadi, R. R. (1998) Large-Scale Simulation of Supersonic Jet Noise, AIAA Paper 98- 
3014, 29th AIAA Fluid Dynamics Conference, Albuquerque, NM. 

6) Choi, D., Barber, T. J, Chiappetta, L. M, and Nishunura, M, (1999) Large Eddy Simulation of High-Reynolds Number Jet Flows, 
AIAA Paper 99-0230, 37th AIAA Aerospace Sciences Meeting, Reno, NV. 


JAMES LOELLBACH 

Research Area: 3D Structured Grid Generation Codes for Turbomacbineiy 

During 1998, James Loellbach performed research tasks relating to numerical analysis of turbomachinery components. The work 
performed in cooperation with Chunhill Hah of NASA Lewis Research Center and Fu-Lin Tsung of ICOMP. 

The work with Fu-Lin Tsung concerned the numerical simulation of flow through the exit volute of a centrifugal compressor for a 
refrigeration system- The simulation was performed using an unstructured Navier-Stokes solver originally developed at NASA Langley 
Research Center (USM3D) and later modified by Tsung, Loellbach and Hah for turbomachinery flows (TUSM). The compressible flow 
through die compressor volute was simulated at a number of different operating conditions to identify the cause of flow separation in the 
exit diffuser. This project is still underway, with the ultimate goal of using the numerical simulations to estimate losses in the volute. 

The work with Chunhill Hah involved grid generation and post processing tasks for numerical simulations of turbomachinery components 
in support of various NASA programs and industry collaborations. These simulations were performed using a structured Navier-Stokes 
solver developed by Hah and grid generation codes developed by Loellbach. One such study examined the effects of blade sweep in a 
transonic axial compressor and yielded information on how to control flow separation near the endwalls. In another study, some simple 
modifications to an existing centrifugal compressor impeller were simulated, revealing ways to improve the efficiency and operating range 
of the compressor. Results of these studies will be published in the near future. 


ANDREW T. NORRIS 

Research Area: Aerothermochemistry 

1. Work Performed 


1.1 Manifold Reaction Schemes. 


Several significant steps were taking in the development of the Intrinsic Low-Dimensional Manifold code. This code takes a full reaction 
mechanism, and automatically simplifies it into a reduced mechanism. 


Work performed consisted of introducing an adaptive tabulation scheme into the code, which provides an automatic way of ensuring a 
more uniform accuracy to the resulting table. This is achieved by placing more grid points in regions of strongly-varying chemical reaction, 
and so achieving better resolution in these regions. 


A second major improvement was the extension of the tables to include super-equilibnum compositions. It is observed that mixing can 
produce chemical impositions that lie "above" the Mly-bumt state, and so the ILDM code was re-written to include these regions. This 
involved utilization of Simplex schemes to establish the limits of composition space, and alteration to the storage methods. 


During this time, the ILDM scheme was used in the NCC code for several different cases. Timing results for the code show significant 
benefits to the performance of NCC compared to using 12 species finite rate methods. 
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No. CPUs 

Seconds per iteration 
12 species ILDM 

Improvement 

Efficiency 

1 

789.2 

147.2 

5.4x 

- 

4 

- 

45.0 

- 

81.8% 

8 

- 

21.9 

- 

84.0% 

16 

52.2 

10.8 

4.8x 

85.2% 

20 

- 

9.2 

- 

80.0% 

24 

- 

7.7 

- 

79.7% 

28 

- 

7.4 

- 

71.0% 

32 

28.5 

6.6 

4.3x 

70.8% 


Implementation of ILDM in the NCC code and a comparison of the results obtained by a 12 step mechanism were reported for the case 
of a piloted methane flame at the 1998 Joint Propulsion Conference. 

1.2 Neural Network Reaction Storage. 

No work was performed on the Neural Network storage of ILDM mechanisms. This part of the project was delayed a year to allow 
Abtech to port their program to the Unix environment. This work is expected to continue in 1999. 

1.3 PDF in NCC. 

Late in 1998, Allison Gas Turbines provided us with a Scalar PDF solver. Currently this solver does not work in a multi-processor 
environment, and this is to be attempted in 1999. 

1.4 Thermodynamic Approach to Reduced Chemical Kinetics. 

Last year, some numerical tools for the investigation of reaction kinetics were developed. Some further work on these was performed 
during the year. Areas investigated included use of minimization techniques to find reaction paths, decomposition of composition space 
into orthogonal basis vectors and comparison with existing kinetic rate expressions. However no significant breakthroughs can yet be 
reported. 

2. Continuing Work 

During the next year, the following projects are proposed: 

1. Extending ILDM method to allow for variable inlet properties. 

2. Validation of the new ILDM scheme in the NCC code. 

3. Validation of the scalar PDF module in the NCC code. 

4. Testing and development of Neural Network approach to storing reduced mechanisms. 

5. Further work on the thermodynamic approach to simplification of chemical kinetic expressions. 
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KRISH RADHAKRISHNAN 

Research Area: Improved Numerics for Reacting Flow Application 

The research work for the period 1/1/98 - 12/31/98 focused on providing analytical and experimental support to the High Speed 
Research, Fast Quiet Engine, Smart Green Engine, Fire Safety, and Rocket Based Combined Cycle ("Trailblazer") programs at the NASA 
Lewis Research Center. The progress made in these different efforts is described below. 

For the High Speed Research, Fast Quiet Engine, and Smart Green Engine Programs, several improvements were made to the NASA 
Lewis kinetics and sensitivity analysis code, LSENS. The ability to handle several new rate coefficient expressions that are being 
increasingly used were included. (The previous version of LSENS could treat only the so-called modified Arrhenius expressions.) The 
debugging of this coding is currendy being done. The automatic calculation of sensitivity coefficients (Le., partial derivatives) of the heat 
release rate with respect to initial condition values and rate coefficient parameters for both static and inviscid, one-dimensional chemically 
reacting flow problems was incorporated. These sensitivity coefficients are required for performing linear stability analysis of combustion- 
acoustic interactions using detailed reaction mechanisms. This new option can be utilhed for the Fast Quiet Engine Program to e xam i n e 
the effects of perturbations in the fuel-air ratio due to fluctuations in die fuel system and/ or the air line. Effects of fluctuation in fuel-flow 
rate due to fluctuations in die fuel line can now be separated from those caused by fluctuations in the air stream. Reeendy, several methods 
have been proposed to automatically reduce a given mechanism, producing a simplified mechanism. One such method, involving the raking 
of a mechanism to discard species and reactions considered unimportant according to criteria provided by the user, was incorporated into 
LSENS. These criteria include for example the contribution of a reaction to the heat release rate and the rate of production of a species 
of interest, such as a pollutant 

Also for the High Speed Engine Research Program, combustion-acoustics interaction data were collected in CE5 for a variety of combustor 
chambers and fuel injectors and for various flow conditions. These results are helpful in identifying operating conditions that result in 
excessive combustion noise or even in resonance, which can damage the combustor. Such studies also enable one to assess the acoustics 
of given combustor/ fuel injector combination and make comparisons among various combustor configurations. 

The following papers were presented in the above areas: 

Yang, Z., Breisacher^C, Oyediran, A., and Radhaknshnan, K., "Combustion Acoustic Instability Analysis of LPP Combustor, Part I: 
Acoustic Vield," AIAA Paper 98-0474, 36th Aerospace Sciences Meeting and Exhibit, Reno, NV, January 12-15, 1998. 

Nguyen, Q.V. and Radhakns hnan , K., "Real-Time Laser Absorption Fuel-to-Air Ratio and Fluctuating Pressure Measuremenents in a Lean 
Premixed Prevaporized Combustor at High Pressure and Temperature," Poster W1E15, Engines Session, 27th International Symposium 
on Combustion, Boulder, CO, August 2-7, 1998. 

Also for die Smart Green Engine Program, the collaboration with Prof. Thomas Hagstrom of the University of New Mexico on developing 
a new flame code that uses very high-order accurate methods was continued. A new adaptive meshing strategy based on minimization (on 
the mesh) of temperature gradients was successfully implemented.. Experiments in one space dimension showed that this new strategy 
reproduces die accuracy of our older mesh generator. The advantage is that our new strategy is applicable in multiple space dimensions. 
High-order time-stepping based on new spectral deferred correction algorithms was also successfully implemented We now have a method 
that is fourth order accurate in space and time, is robusdy stable for complex reaction and diffusion mec han i sm s, and allows the use of 
Steady simph&ed diffusion operators in the implicit solver. We are very hopeful that these methods will allow us to use simple dimensional 
splittings in higher space dimensions while still providing arbitrary order of accuracy. We have completed a thorough testing of die 
algorithm via computations of propagating plane flames and comparisons with results obtained with our plane flame code— fourth order 
convergence under refinement was verified. We have also established the accuracy of our code for plane steady flames by comparing the 
calculated laminar flame speeds with experimental data for a variety of fuels and initial conditions. 

The following papers were presented: 

Hagstrom, T., K-, and Zhou, R, "Computation of Steady and Unsteady Laminar Flames: Theory," AIAA Paper 98-3246, 

34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Cleveland, OH, July 13-15, 1998. 

Hwang, S.M., Radhaknshnan, K., Nguyen, Q.V., and Rabinowitz, M.J., "Is there an S04 Intermediate in S02 Oxidation?" Poster W4B16, 
Reaction Kinetics Session, 27th International Symposium on Combustion, Boulder, CO, August 2-7, 1998. 

Hagstrom, T., Radhaknshnan, K-, and Zhou, R., "Efficient and Reliable Computation of Steady and Unsteady Flames," Poster W5D 1 1 , 
Premixed Models and Simulations Session, 27th International Symposium on Combustion, Boulder, CO, August 2-7, 1998. 
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For the fire safety program, die effort during this period consisted of literature review of fuel properties, flammability, and testing. Also, 
along with Bryan Palaszewski, the needs for laboratory testing were defined and an in-house experimental program for measuring 
flammability limits and minimum ignition energy of aviation fuels was planned Discussions were also held with Boeing personnel in order 
to define a model problem for die CFD (computational fluid dynamics) simulation of aircraft fuel tanks. 

For the trailblazer program, the equilibrium routines built into LSENS were rewritten and adapted for use with a 1-D CFD code, 
developed by Dr. Shaye Yungs ter for analyzing the "independent ramjet stream" (ERS) cycle, in which the rocket and ramjet streams do 
not mix. Also for the hypersonics program the collaboration with Dr. Yungs ter on Detonation Wave Modeling continued. The goals of 
this work are to continue establishing validity of our two-dimensional, time-accurate CFD code by application to standing oblique 
detonation waves for which experimental data are available and to improve understanding of detonation wave physics. In this study, we 
investigated the flow of a hydrogen-air mixture over two different projectile configurations in an expansion tube and comparisons made 
with OH PLIF data. These CFD studies complement experimental work by providing detailed information about combustion initiation, 
flow structure and flow establishment time, a critical parameter in pulse facilities, in which the available test time may be too short to 
establish fully the reacting flowfield. Such calculations are also appropriate benchmark cases for evaluating numerical models for flow and 
combustion chemistry. 

The following papers were published/presentcd: 

Yungster, S., Radhakrishnan, K., and Rabmowitz, M.J., "Reacting Flow Establishment in Ram Accelerators: A Numerical Study," Journal 
of Propulsion and Power, VoL 14, No. 1, pp. 10-17, January-February 1998. 

Yungster, S. and Radhakrishnan, K., "Simulation of Unsteady Hypersonic Combustion," Poster W4C 13, Energetic Materials and Supersonic 
Combustion Session, 27th International Symposium on Combustion, Boulder, CO, August 2-7, 1998. 

AAMIR SHABBIR 

Research Ana: Turbulence Modelling 

The arithmetical average surface roughness height for a typical compressor blade ranges from 20 to 80 micro-inches. The Reynolds 
number based on this length scale could range from a transitional to a completely rough regime, especially in the last blade rows of a 
multistage turbomachine. This means that for such cases the skin friction obtained by assuming a hydraulically smooth surface will be 
deficient. 

The effect of surface roughness on skin faction is incorporated in a wall function approach which uses Spalding's formula. The skin 
friction calculation is an iterative process since the unknown appears on both sides of the equation in an implicit way. A guess of friction 
velocity is required before a new value of skin friction can be calculated. The procedure is repeated until the new value of skin friction 
satisfies some convergence criteria. This procedure can diverge in certain situations (say a bad initial guess for starting a simulation). Thus 
from the CFD point of view an explicit procedure is preferable. 

Therefore, a methodology is formulated for explicit calculation of skin friction through the use of an additional variable - the ratio of 
the surface roughness height to the distance from the solid surface. For a given grid this information is readily available in a CFD code. 
This methodology is incorporated into a multistage turbomachinery code. 

SCOTT SHERER 

Research Area: Simulations of Electromagnetic Phenomenon 

The work that Dr. Sherer performed in the Computational Sciences Branch, Air Force Research Laboratory, Wright-Patterson AFB, during 
this period falls into three main categories. These categories are interrelated and a brief summary of each is given below. 

1. CHARGE Validation: Several validation cases were studied using “CHARGE”, the Finite-Volume Time-Domain Maxwell 

equation solver developed by the branch as part of the DoD High Performance Computer Modernization Office’s (HPCMO) 
Common High Performance Computing Software Support Initiative (CHSSI). The three-fold purpose of these studies was to 
validate the flux-vector splitting/Runge-Kutta algorithm used in the code against a wide variety of benchmark geometries, to 
establish a centralized reference and baseline performance data for future work, and to directly compare the accuracy and 
efficiency of CHARGE to predict radar cross-section (RCS) signatures against a popular and efficient frequency-domain 
computer program. Several geometries, including spheres, flat plates of various platforms, cubes, ogive cylinders and finned 
missiles were examined. Grid refinement studies and comparisons of various topologies for some of these geometries were also 
performed. This work will be presented in a paper Dr. Sherer is co-authoring at the AIAA Aerospace Sciences Meeting. 
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2. Further Development and Enhancement of CHARGE: Modifications were made to CHARGE to enhance its capabilities. 

Foremost among these enhancements was adding to the ability to model a Gaussian pulse incident wave scattered by a body. 
By using a Gaussian pulse, which contains information over the frequency spectrum, RCS profiles at various frequencies can 
be extracted from a single simulation with greatly reduced resources expended. This capability was successfully applied to 
calculate RCS profiles for multiple frequencies for a sphere, which compared very well with single frequency predictions as well 
as theoretical data. 

3 Prediction of RCS for Full Bl-B Aircraft: Dr. Sherer also performed extensive grid generation work in support of the 

Computational Sciences Branch’s DoD HPCMO Challenge Project entided “B-1B Radar Cross-Section Prediction”. This gnd 
generation involved siting die existing grid, which consisted of over 700 computational blocks, and consolidating this number 
down to a more f^nageahle size. Concurrent, various modifications to die grid itself were made to improve the overall quality 
of die grid as well as increase the minimum allowable time step, which would result in shorter execution times. These various 
n ,^tifi^rinnc yielded mixed results. The final report on this project will be presented to the HPCMO. 

SHYUE-HORNG SHIH 

Research Area: Direct Prediction of Supersonic Jet Noise 

In supersonic jets, both large-scale coherent structures and small-scale turbulent eddies contribute to the radiated noise in the far field. 
The existing Large Eddy Simulation code was applied to the test case of a Mach 1.4 heated round jet in a 400 ft/s uniform free stream. 
Unlike die previous applications of the LES code, the large-scale coherent structures are not the dominant noise generation mechanism 
in this jet The calculation of this jet will aid in determining the ability of the algorithm to capture turbulent noise generation as well as 
instability wave noise generation. The simulated jet was excited with multiple frequency disturbances at the jet plume inflow. The far fiel 
noise was by n*mg Kirchhoff s method with the input information coming from the near field LES results. The computed mean 

flow profiles and noise directivity in the far field compared favorably with the experimental data. However, the noise spectra show that 
the r°-" f calculation was unable to capture die high frequency regime of noise radiation. To resolve this issue, further investigation is 
underway. One approach would be to increase the mesh resolution and another is to use an advanced, high-order accuracy, numencal 
scheme for the LES code. 

To the c«mp»wrinnal costs for die LES method, an alternative approach has also been explored. Based on the LES technique, the 

fluid is split into three kinds of motion: a time-averaged motion, a large scale wavelike structure, and small-scale, random, 

taAsaleace Non-linear disturbance equations, which are derived from the Navier-Stokes equations, are solved with the prescribed mean 
flow The lame s cale structures are computed direedy while the effects of small-scale turbulence are modeled. The full three dimensio 
calculation shows that the development of the disturbances are in good agreement with experimental observations, while the computer 
CPU time required is 50% of that of the full LES simulation. It should be noted that in this approach the mean flow has to be specified 
and the high frequency resolution issue in the full LES approach is also encountered. 

TSAN-HSING SHIH 

R / search Area : Turbulence Modelling 

The /-»H^ nonlinear turbulence model has been refined and tested for swilling flows which have critical importance in die design of aircraft 
engines. 

1 The model was upgraded to be tensorially invariant so that there is no need to define and calculate the wall distance for each point 
in the computational domain. For generated numberical grids, this will save hours of CPU time in each preprocessing . 

2. The model was evaluated for other benchmark problems and separated flows. 

i) In backward facing -step flows, the model was able to simulate die fine flow structure in the comer region. A senes of 
counter-rotating eddies was predicted and verified by comparison with experimental data. 

ii) In Roback-Johnson combustor flows, the center recirculation bubble and the wall separations were accurately predicted. 
The swirling velocity in the central region of the combustor was also predicted well, as verified by comparison with existing 
experimental data. 

iii) The model also worked well in other traditional benchmark flows: flat plate boundary layers, channel flows, mixing layers 
and jets. 
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3. A general scalar equation was successfully included in a two-dimensional Naviex-Stokes code. This will enable die evaluation of a 
turbulent heat transfer model for complex flows. This will become a research project for 1999. 

The research results were partially reported in the following AIAA papers: 

1. AIAA 98-2551: “A k-€ Model for Wall-Bounded Shear Hows.” 

2. AIAA 98-3983: “A Nonlinear k-E Model for Turbulent Shear Hows.” 

3. AIAA 98-3243: “Turbulence Model Developments at ICOMP.” 

ERLENDUR STEINTHORSSON 

Research Area: Code Development for Flows tn Complex Geometries such as Turbine Blade Coolant Passages 

During the last year, the TRAF3D.MB code (a multi-block and multi-grid flow solver for internal flows) was modified to allow parallel 
execution on shared-memory parallel computers (SMP computers). This was done through the use of compiler directives that enable the 
compiler to parallelize the execution of loops that include subroutine calls. With this approach, each CPU of the parallel computer handles 
one grid block at a time, performing all the compute-intensive operations on that block. Synchronization occurs after all blocks have been 
advanced one stage in the Runge-Kutta time stepping scheme that is used in die code to advance the solution in time to steady state. Since 
the SMP model was used, no special treatment was required for exchanging data between blocks that are handled by the different 
processors. The code has been tested on CRAY -CSX) and CRAY-J90 computers at NAS and has been found to work well. No significant 
overhead due to parallel computations is experienced since the CRAY is operated as a time-shared system. Thus, during parts of the 
program that are not executed in parallel, die idle processors are utilized by other users. The code was also ported to SGI 0rigin2000 
systems at NAS. The Origin 2000 computer systems differ form the CRAY systems in that memory is distributed on the processor nodes. 
Thus, it is not a true SMP computer in hardware. However, the operating system and compilers allow die user to program the computer 
as if it were an SMP. In this case it was found that for small problems, the parallelization strategy that worked on die CRAY systems also 
worked for the SGI systems. For larger problems, the execution fails on multiple CPUs. Even with expert analysis from NAS, Even with 
expert analysis from NAS, this issue has thus far not been resolved and is unexplained at the present time. 

During die last year, a detailed report was written that describes the design of the TRAF3D.MB code as well as the governing equations 
and discretization used in the code. The report will be published by NASA Lewis Research Center (now J. Glenn Research Center at Lewis 
Field). The TRAF3D.MB code was officially made available to the gas-turbine engine community at an IHPTET meeting in September 
1998. 


GERALD TRUMMER 

Area : System Administration 


Maintained and upgraded the operating system software for the AFRL/VAA workstation cluster. This included installing new software, 
software patches and editing and configuring files to bring new machines into the network cluster and making the machines secure. 
Performed other system administration duties such as archiving user accounts, adding/removing users, backing up file systems, installing 
new hardware, etc. 

Installed, maintained and upgraded visualization plotting packages used by the VAA division. This included new packages requiring 
configuration of license servers and script files and editing existing software at users' request to add or improve features. 

Made scientific visualization movie of Drs. Visbal and Gordniers* “Impingement of a Streamwise Vortex on a Plate” to be submitted to 
American Physical Society's Gallery of Huid Motion. Made sci-vis movie of Dr. Ge-Cheng's Inlet How for presentation at a conference. 
Made other sci-vis movies for the VAA division for presentation at various conferences. Made multi-media web page of Dr. Rizzetta's 
'Numerical Investigation of Synthetic Jet Howfields' AIAA paper. http://www.cfd.wpafb.af.mil/Pubs/vaac/rizzetta/synjet. Maintained 
other web pages and the web server. 

FU-LIN TSUNG 

Research Area: Development of 3D Structured! Unstructured Hybrid Navier-Stokes Solver for Turbomacbineiy 

The goal for the present research is to develqp, validate, and apply a 3-D, Navier-Stokes, unstructured-grid solver for turbomachinery 
applications. For 1998, the primary development phase of the solver is complete for vector machines. The solver has been applied to an 
unsteady internal flow passage to study the effects of cavity seals and the solution has been distributed to interested industrial researchers. 
The solver is also being used to study the performance of a complex volute geometry. The flexibility of this unstructured system was 
demonstrated in the volute case by its ability to rapidly return preliminary results of the complex flow field, from grid generation to solution 
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some 
a more accurate 


Within a week For the preliminary analysis, a coarse grid of 182,000 cells was used, which translates to 35 MW of memory and requires 

,5-30 C90ho«= 

4 processors, the gain drops off quickly and diminishes the value of using additional processing resource. 

ELI TURKEL 

Research Ana- PnconditioningApplud to Turbomacbtneiy Plow Simulation 

We worked on the development of preconditioning matrices for compressible low speed flows. Work has been for 

applications to external flows. As part of the work at ICOMP, we considered applications GRC ^S^were 

narticular the preconditioning was incorporated in two different internal flow/ turbomachinery codes at NASA GRC. Applications w 
Skfh« involve a threHLensional stator and stator-rotor configurations. These results demonstrated the utility of ^ 

^ iniTof convergence acceleration and in terms of accuracy for low speed flows. It was found that it was especially important to 
account for the inflow boundary conditions without resorting to characteristic variables. 

In addition work was done on extending the SLIP/CUSP scheme for an artificial viscosity to internal codes. This also required 
changes in the original external flow version of these algorithms. This was also combined with the precondmong to yield a more ac 
treatment of boundary and shear layers. 

A third project was die construction and analysis of a compact implicit extension of the MacCormackschemeand 

variants. This allowed the development of an efficient code for three dimensional acoustics with fourth and sixth order accuracy. This 

joint work with Dr. Ray Hixon from ICOMP. 

MICHAEL WHITE 

Research Area: Electromagnetics and Fluid Dynamics 
development of the capabilities of the code. 

not consistent with the original PDEs and the nature of this inconsistency was investigated. A paper on the topic is in progres 
Virginia 121 . 

Dmg db, «, . .ho* ™ ^ » Udfcing . 

u G«»™h — *» -d ■**« » “ 

functionality of program and of the user interface to the program. 


References 


[1] 

[ 2 ] 


White, M. D. and Kollmann, W. "Consistency of Filters for Stabilizing Numerical Solutions of PDEs", (1999) in 
preparation. 

White, M. D. and Kollmann, W. "Hybrid High-order Compact Finite-difference/Spectral Scheme for Simulanon of Round 
Jet Rows", AIAA Paper 99-3258 
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SHAYE YUNGSTER 

Research Area: High Spetd Combustion and Detonation Waves 
The research work for the 1998 year focused on two different areas described below 
1) CFD Analysis of Rocket-Based Combined Cycle (RBCQ propulsion systems at low speed 

The NASA Lewis Research Center is currently developing a reusable, single-stage to orbit launch vehicle known as "Trailblazer" that is 
based on a Rocket-Based Combined-Cycle (RBCC) propulsion system- This vehicle will operate in four modes from lift-off to orbit 1) 
ejector-ramjet, 2) ramjet, 3) scramjet, and 4) all-rocket 

An initial study of the all-rocket mode (described in the 1997 ICOMP report) was completed during the first quarter of 1998, and the 
results were s ummariz ed in the following publications: 

C. J. Steffen, C.J., Smith, TD., Yungs ter, S. and Keller, D .J., "Rocket Based Combined-Cycle Analysis Using NPARC", AIAA paper 98- 
0954, January 1998. 

C.J., Smith, TX>., C. J. Steffen, Yungs ter, S. and Keller, D J., "Performance of an Axisymmetdc Rocket Based Combined Cycle Engine 
During Rocket Only Operation Using Linear Regression", NASA TM- 1998-206632, March 1998. 

C J-, Smith, T.D., C. J. Steffen, Yungster, S. and Keller, D .J., "Analysis of Rocket Based Combined Cycle Engine During Rocket Only 
Operation", NASA TM-1 998-206639/REV 1 October 1998. 

The subsequent research work focused on the low speed operation mode, that is the ejector-ramjet, which is described below. This mode 
of operation typically covers the Mach number range from takeoff to approximately Mach 3. 

In the conventional ejector-ramjet operation mode, a fuel-rich rocket exhaust is mixed and burned with air captured by the inlet The rocket 
provides all of the fuel needed for combustion with the entrained air. The internal flowpath is designed to produce thermal choking where 
mixing is complete. The products of combustion subsequently expand through the vehicle nozzle. The main disadvantage of this concept 
is the relatively long duct required to achieve complete mixing of the air and rocket streams. In order to overcome this difficulty, a variation 
of the conventional ejector-ramjet was proposed in which the requirement for complete mixing of die two streams is removed. In this new 
"Independent Ramjet Stream" (IRS) cycle, the airstream is fueled independently using the ramjet and scramjet mode fuel injectors. This 
can be accomplished upstream, since the air stagnation temperatures during this mode are not high enough to cause autoigmtion. The 
rocket will serve as a pilot for the fueled airstream. A single rocket can therefore be used without a long mixing duct, and since the rocket 
is not the fuel source for the airstream the rocket oxidizer/ fuel (O/F) ratio can be scheduled for best system performance. In addition, 
the fuel injectors provide the to control the location of the thermal choke by adjusting the amount and penetration distance of the 

fuel injected into die airstream. 

Calculations have been performed to analyze both of these ejector-ramjet concepts, however most of the effort has been recently directed 
towards the IRS cycle analysis. Two approaches have been taken; 1) development of a quasi-one dimensional CFD code for fast analysis 
of the performance of the IRS cycle, and 2) application of an axisymmetnc CFD code to provide a more detailed analysis of the flow and 
combustion processes. 

Quasi-one dimensional CFD calculations of the IRS cycle 

The objective of this work is to develop an analysis tool for the IRS cycle, investigate the effects of airstream fuel-air ratio, heat release 
distribution, and rocket flow conditions on system performance over the operating flight range of this mode of operation. An additional 
objective is to generate performance maps for trajectory optimization. 

Approach 

The ramjet and rocket streams are solved simultaneously using a TVD MacCormack time-marching scheme. A quasi-one dimensional 
approximation is used to model both streams. Combustion in the ramjet stream is modeled by a prescribed distribution of hydrogen fuel 
along the combustor duct Equilibrium chemistry is used to model the combustion process utilizing the LSENS kinetics code of 
Radhaknshnan. The two streams are coupled using an auxiliary equation that adjusts die cross-sectional area of the streams such that the 
static pressures are equal at all points, and the sum of the areas equal to the available duct area- 


16 




RESEARCH IN PROGRESS 


Accomplishments 

This code has been applied for generating performance maps for a specific trailblazer configuration covering the entire range : of 1 Mach 
numbers (0 - 3) and rocket chamber pressures (100-1500 psi> The effects of varying the airstteam fud-air raUo was investigated and die 
flowpath was adjusted for optimum performance. A summary of this work will be presented at die 35th Joint Propulsion Conference. 

A ti symmetric CFD calculations 

In this work, we are analyzing an axisymmetnc configuration that closely models the trailblazer vehicle operating on the IRS cycle or the 
conventional ejector-ramjet mode. Most of the recent work has focused on the IRS cycle, and die discussion below concentrates on dus 
propulsion mode. The goal is to understand the flow and combustion physics, study the effects of airstream ftiel-air ratio, mixture 
distribution, geometric configuration, and rocket chamber pressure and O/F ratio on engine performance 

Approach 

The analysis of the IRS cycle is carried out using a specialized reacting flow code developed in-house. This code has been ^ed m the past 
to accurately model combustion phenomena in high-speed propulsion applications. It sokes the Navier-Stokes equauons wi* 
chemistry and real gas effects using an implicit, total variation diminishing (TVD) algonthm. It includes a generalized detailed . 

capability, various options for turbulence models, and steady-state or time accurate marching algonthms. The hydrogen-air combustion 
is modeled with a 7-species, 8-step reaction mechanism. 

Accomplishments 

Calculations have been performed to examine the effects of airstream ftiel-air ratio, mixture distribution, and rocket cha^ pressure on 
flame propagation and stability over the operation range of this propulsion mode (from see-level static conditions to 

SpJL Impulse, heat release distribution and other performance parameters have been obtained for vanous ejector-rmxer configmanons. 
Ite^alrito demonstrate thermal choking and control of thermal choke location by varying the properties of the fueled airstream. 

This detailed analysis of the IRS mode of operation is required in order to determine system performance develop 

use in engine cydelodes, and support experimental efforts. In particular, planned experimental testing ofRBCC configurations NASA 
GlenmResearch Center will rdy on the computational analysis for determining the conditions and configurations to be run in the 
experiments, and for post-test analysis of the data. 

2) Detonation Wave Modeling 

We are continuing work on detonation wave modeling. The goal of this work, which is carried out in collaboration with K. Radhaknshnan 
rito both continue establishing validity of the CFD code we have developed, by application to standing oblique detonauon waves for which 
experimental data has been recendy reported, and improve understanding of detonation wave physics. In this study, we investigated the 
fl^^hydrogen-air mixture over twVdifferent projectile configurations in an expansion tube. The computed solutions will be com f^' d 
iith JperSToHPLIF data. These computational studies complement experimental work by providing detailed ' f*™ 
combustionminaaon flow structure and flow establishment time, a critical parameter in pulse facilities, in which the available test time 
may be too short to establish ftilly the reacting flowfield. Such calculations are also appropriate benchmark cases for evaluatmg numencal 
models for flow and combustion chemistry. Results of this work were reported in two papers: 

Yungster S., K., and Rabinowitz, M.J., "Reacting How Establishment in Ram Accelerators: A Numerical Study," Journal 

of Propulsion and Power, Vol. 14, No. l,.pp. 10-17, January-February 1998. 

Yungster, S. and Radhakrishnan, K., "Simulation of Unsteady Hypersonic Co^ustion," Poster W4C13 Energetic Materials and Supersonic 
Session, 27th International Symposium on Combustion, Boulder, CO, August 2-7, 1998. 

Additional results will be presented at the 35th Joint Propulsion Conference. 

GE-CHENG^ZHA^ ^ Navier-Stokes Time Accurate Solutions Using Mubtpartitiomng Parallel Computation Methodology 

A narallel CFD code solving 3D time accurate Navier-Stokes equations with multipartionmg parallel Methodology *being developed m 
i^nhnwT' oSe University within the Air Vehicle Directorate, at Wright Patterson Air Force Base. The advantage of he 
m ^upv6ot>mg parallel method is hat die domain decompostion will not introduce domain boundaries for the implicit operators. g 
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structure data communication is employed so that the im p licit time accurate method can be implemented for multi-processors with the 

same accuracy as for the single processor. No sub-iteration is needed at the domain boundaries. * 

The code has been validated for some typical unsteady flows, which include Coutte Flow, flow passing a cylinder. The code now is being 
employed for a large scale time accurate wall jet transient flow computation. The preliminary results are promising. The mesh has been 
refined to capture more details of the flow field. Hie mesh refinement computation is in progress and would be diffic ult to successfully 
implement without the parallel computation techniques used. A modified version of the code with more efficient inversion of the 
diagonalized block matrix is currently being tested. 

MEI ZHUANG 

Rtsearcb Ana : 

Acoustic problems arise in many engineering applications such as jet noise, the duct acoustic, etc. As we can see, these aeroacoustic 
problems are time dependent, the noise frequency is generally very high, and the noise spectrum is fairly wide. These require that the 
numerical schemes be almost free of dissipation and dispersion in order to ensure that waves travel with an accurate speed, phase, and 
amplitude. 

Several benchmark aeroacoustics problems are studied using a new high-resolution and multi-dimensional numerical method developed 
at NASA Lewis 11,21 . The potential applications of such a scheme in practical engineering problems will be investigated. 

References 

1. S. C. Chang, ‘The Method of Space-Time Conservation Element and Solution Element Method - A New Approach for Solving 
the Navier-Stokes and Euler Equations/’ J. Comput Phys v, 119, 1995, pp. 295-324. 

2. S. C. Chang, X. Y. Wang, and C. Y. Chow, ‘The Space-Time Conservation Element and Solution Element Method - A new 
High-Re solution and Genuinely Multidimensional Paradigm for Solving Conservation Laws - I. The Two-Dimensional Time 
Marching Schemes/’ submitted for publication in J. Comput. Phys., 1998. 
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1998 REPORTS AND ABSTRACTS 


Keith Theo G , Balog, Karen, and Povineffi, Louis A., Editors: “Institute for Computational Mechanics in Propulsion 
(ICOMP), Twelfth Annual Report - 1997”, TCOMP Report 98-1 , NASA TM 208397, June, 1998, 44 pages. 

The Ins titute for Computational Mechanics in Propulsion (ICOMP) is operated by the Ohio Aerospace Institute 
(OAI) and funded under a cooperative agreement by the NASA Lewis Research Center in Cleveland, Ohio. The 
purpose of ICOMP is to develop techniques to improve problem-solving capabilities in all aspects of computational 
mechanics related to propulsion. This report describes the activities and accomplishments during 1997, the Institute’s 
twelfth year of operation. 

Choo, Yung K., and Slater, John W. (NASA Lewis Research Center); Henderson, Todd L. (Raytheon Systems 
Company); Bidwell, Colin S., and Braun, Donald C. (NASA Lewis Research Center); Chung, Joongkee (ICOMP): 
“User Manual for Beta Version of TURBO-GRD: 2A Software System for Interactive Two-Dimensional Boundary /Field 
Grid Generation, Modification, and Refinement”, ICOMP Report 98,2, October, 1998, 67 pages. 

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation, modification, 
and refinement Its features allow users to expliddy control grid quality locally and globally. The grid control can 
be achieved interactively by using control points that the user picks and moves on the workstation monitor or by 
direct stre tching and refining. The techniques used in the code are the control point form of algebraic gnd 
generation, a damped cubic spline for edge meshing, and parametric mapping between physical and computational 
domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry 
manipulation. Internal block boundanes are constructed and shaped by using a Bezier curve. Because TURBO- 
GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the backgroun 
of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can 
be called an interactive solution-adaptive grid generation. 


Hagstrom, T. (The University of New Mexico and ICOMP); and Hariharan, S. I. ^ ^mversity ofARron): “A 
Formulation of Asymptotic and Exact Boundary Conditions Using Local Operators”, ICOMP Report 9& U, NAS 
207935, June 1998, 22 pages. 

In this paper we describe a systematic approach for constructing asymptotic boundary conditions for isotropic wave- 
like equations using local operators. The conditions take a recursive form with increasing order of accuracy^ In three 
dimensions the recursion terminates and the resulting conditions are exact for solutions which are descnbed by finite 
combinations of angular spherical harmonics. First, we develop the expansion for the two-dimensional wave 
equation and construct a sequence of easily implementable boundary conditions. We show that m three dimensions 
and analogous conditions are again easily implementable in addition to being exact. Also, we provide extensions 
of these ideas to hyperbolic systems. Namely, Maxwell’s equations for TM waves are used to demonstrate the 
construction. Finally, we provide numerical examples to demonstrate the effectiveness of these conditions or a 
model problem governed by the wave equation. 

Ameri AJL (AYT Corporation); Steinthorsson, E. (ICOMP); and Rigby, David L. (NYMA, foe.): “Effeccof Tip 
cILg « fL Transfer and Stage Efficacy in A=al Turb»e,-. KWR.pW g K NASA CR 

208514, August, 1998, 18 pages. 

Calculations were performed to assess the effect of the tip leakage flow on the rate of heat transfer to bMe by e 
tip and casing. The effect on exit angle and efficiency was also examined. Passage geometries with and without 
casing recess were considered. The geometry and the flow conditions of the GE-E first stage ^me ^ich 
represents a modem gas turbine blade were used for analysis. Clearance heights of O/o 1 /o, 1.5 /o and 3 A of the 
passage height were considered. For the two largest clearance heights considered, different recess dep ^sw ere 
studkd TbL was an increase in the thermal load on all the heat transfer surfaces considered due to enlargement 
of the clearance gap. Introduction of recessed casing resulted in a drop m the rate of heat transfer on the pressure 
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side but the picture on the suction side was found to be more complex for the s mall er tip clearance height 
considered. For the larger tip clearance height the effect of casing recess was an orderly reduction in the suction 
side heat transfer as the casing recess height was increased. There was a marked reduction of heat load and peak 
values on the blade tip upon introduction of casing recess, however only a small reduction was observed on the 
casing itself. It was reconfirmed that there is a linear relationship between the efficiency and the tip gap height It 
was also observed that the recess casing has a small effect on the efficiency but can have a mnH^ting effect on the 
flow underturning at smaller tip clearances. 

Hixon, R. (ICOMP): “Nonlinear Comparison of High-Order and Optimized Finite-Difference Schemes”, ICOMP Renort 
98-5. NASA CR 208670, October, 1998, 23 pages. 

The effect of reducing the formal order of accuracy of a finite-difference scheme in order to optimize its high- 
frequency performance is investigated using the 1-D nonlinear unsteady invisdd Burger’s equation. It is found that 
the benefits of optimization do carry over into nonlinear applications. Both explicit and compact schemes are 
compared to Tam and Webb s explicit 7-point Dispersion Relation Preserving s chem e as well as a Spectral-like 
compact scheme derived following Lele’s work. Results are given for the absolute and L2 errors as a function of 
time. 

Hixon, R., Shih, S.-H. (ICOMP); Mankbadi, R. R. (ICOMP and Cairo University): “Evaluation of Boundary 
Conditions for the Gust-Cascade Problem”, ICOMP Report 98-6. NASA CR 208671 , November, 1 998, 34 pages. 

Using a high-order accuracy finite-difference time-domain algorithm, the acoustic scattering from a flat-plate cascade 
is computed Keeping the gnd and time step fixed, the effect of four different boundary conditions on the accuracy 
and stability of the computed solution is compared. 

Hixon, R. (ICOMP); Tutkel, E. (ICOMP and Tel Aviv University): “High-Accuracy Compact MacCormack-Type 
Schemes for Computational Aeroacoustics”, ICOMP Report 98-7. NASA CR 208672, October, 1998, 37 pages. 

Using MacCormack-type methods, a new class of highly accurate compact MacCormack-type schemes is derived 
which does not require a tridiagonal matrix inversion to obtain the spatial derivatives. Two examples are shown, 
and results of these schemes for three linear and nonlinear CAA Benchmark Problems are presented. 
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AIAA-98-3240 

Nonlinear Comparison of High-Order and Optimised Numberical Schemes 

R. Hixon, Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH 

The effect of reducing the formal order of accuracy of a finite-difference scheme in order to optimize its high- 
frequency perfo rman ce is investigated using the 1-D nonlinear unsteady inviscid Burgers’ equation. It is found that 
the benefits of optimization do carry over into nonlinear applications. Both explicit and compact schemes are 
compared to Tam and Webb’s explicit 7-point Dispersion Relation Preserving scheme as well as a Spectral-like 
compact scheme derived following Lele’s work- Results are given for the absolute and L2 errors as a function of time. 


AIAA-98-3241 

Time Domain Solution of the Airfoil Gust Problem Using High Order Compact Scheme 

R. Hixon and S.-H. Shih, Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH; R. 
Mankbadi, Cairo University, Cairo, Egypt; J. Scott, NASA Lewis, Cleveland, OH 

A finite difference time domain solution of the airfoil gust problem is obtained using a sixth-order compact scheme 
to compute the spatial differences. For computational efficiency, the equations are cast in chain-rule curvilinear form, 
and a multiblock structured solver is used. Preliminary results are shown for symmetric and cambered Joukowski 
airfoils. These results are compared to solutions obtained by the GUST3D frequency-domain solver. 


AIAA-98-3242 

Effect of Airfoil Ice Shape Smoothing on the Aerodynamic Performances 

J. Chung, Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH; A. Reehorst, Y. Choo, 
and M. Potapczuk, NASA Lewis, Cleveland, OH 

Complicated and irregular ice shapes accreted on an airfoil of turbo-prop aircraft were modelled with various levels 
of smoothness to study its effect on aerodynamic properies Curves with three different smoothing levels were 
constructed by assigning three different percentages of control points (CPs) for selected four representative ice shapes 
obtained by experiment in Icing Research Tunnel at NASA Lewis. Numerical result showed that as the modelled 
surface gets smoother by the use of smaller number of CPs, all the ice shapes indicated a trend which can be 
characterized as increase in lift and decrease in drag. The level of smoothing obtained by using 50% of CPs was 
found to be acceptable in most cases, returning very close results compared to that of 100% CP level. Among the 
ice shapes examined, a prominent redged-ice formation was found to cause the largest aerodynamic penalty. 


AIAA-98-3243 

Turbulence Model Development for Propulsion System at ICOMP 

T.-H. Shih, Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH 

Some recent model developments at ICOMP, NASA Lewis Research Center are briefly described. The discussion 
covers a wide range of models from Reynolds-stress algebraic equation models to Reynolds -stress transport equation 
models. Some validations and applications performed by ICOMP researchers are presented for each models described 
in this paper. 


AIAA-98-3244 

Modeling of Shock Wave/ Turbulent Boundary-Layer Interaction 

W. Liou, Western Michigan University, Kalamazoo, MI, T.-H. Shih, Institaute for Computational Mechanics in 
Propulsion, Cleveland, OH 
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AIAA-98-3983 

A Non-Linear k-e Model for Turbulent Shear Flows 

T.-H. Shih, Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH; K.-H. Chen, The 
University of Toledo and NASA Lewis Research Center, Cleveland OH; N.-S. Liu, NASA Lewis Research 
Center, Cleveland, OH 

The model introduced in this paper is an extension of the linear k-e eddy viscosity model (Shih et al, 1995) to a cubic 
stress-strain relation (Shih et al, 1997) and can be applied down to die solid walL The model is shown to perform 
very well for many benchmark shear flows including wall-bounded flows (e.g., flat plate boundary layers, channel 
flows). It is also shown that the effect of swirling on turbulent flows can be effectively modeled by the present model 
Another feature of this model is that the model does not contain wall distance or any wall related parameters, so that 
the model is tensodai invariant and also Galilean invariant. This feature of the model is very useful to various 
engineering applications for flows with complex geometries and to CFD codes with unstructured grids. 
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Computation of Steady and Unsteady Laminar Plane Flames: Tbeoiy 

K. Rad h a k ris h an an , Institute for Computational Mechanics in Propulsion (ICOMP), Cleveland, OH; T. 

Hagstrom, University of New Mexico, Albuquerque, NM 

In this paper we describe the numbencal analysis underlying our efforts to develop an accurate and reliable code for 
si mu lati n g flame propagation using complex physical and chemical models. We discuss our spatial and temporal 
discretization schemes, which in our current implementations range in order from two to six. In space we use 
staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion 
operators while m ai ntainin g ellipbcity. Our temporal discretization is based on the use of preconditioning to produce 
a highly efficient lineady implicit method with good stability properties. High order for time accurate simulations is 
obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for 
computing stationary flames. The primary issue here is the automatic generation of initial approximations for the 
application of Newton’s method. We use a novel time-stepping procedure, which allows the dynamic updating of 
the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, 
primarily for die stationary flame problem. These illustrate the reliability of our techniques, and the dependence of 
the results on various code parameters. 
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Time Accurate 3D Navter-S takes Analysis of a Subsonic Axial Turbine Stage 

H. Bcnetschik, AEA Technology, Otterfing, Germany; T.-W. Volmar, B. Brouillet, and H.-E. Gallus, Institute fur 
Strahlantriebe und Turboarbeitsmaschincn, Aachen, Germany 
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